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abstract: In this article, we present the properties of the C-parameter which is one of event shape variables. 
We obtain the coupling constants both in the perturbative and in the non-perturbative part of the QCD theory. 
To achieve this we fit the dispersive model as well as the shape function model with our data. Our results are 
consistent with the QCD predictions. We explain more features of our results in the main text. 
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I. INTRODUCTION 


Event shape variables in e'^e~ annihilation provide an ideal 
testing ground to study Quantum Chromo-Dynamics (QCD) 
and have been measured and studied extensively in the last 
three decades. In particular, event shape variables are interest¬ 
ing for studying the interplay between perturbative and non- 
perturbative dynamics m One of the most common and suc¬ 
cessful ways of testing QCD has been by investigating the dis¬ 
tribution of event shapes in e^e~ hadrons, which have been 
measured accurately over a range of centre-of-mass energies, 
and provide a useful way of evaluating the strong coupling 
constant a^. The main obstruction to obtaining an accurate 
value of as from distributions is not due to a lack of precise 
data but to dominant errors in the theoretical calculation of the 
distributions. In particular, there are non-perturbative effects 
that cannot yet be calculated from first principles but cause 
power-suppressed corrections that can be signicant at experi¬ 
mentally accessible energy scales O. 

In this article, we show the cross section for C- parameter as 
an event shape observable. Also we peruse both perturbative 
and non-perturbative theory for calculating coupling constants 
using different models .It is mentioned that we have already 
done some analyses on some event shape observables in our 
previous publications (SJIH. 

The outline of the paper is as the following: In sect.2 we 
define and review the C-parameter and show the cross sec¬ 
tion distribution for different energies. In sect.3, we present 
the calculations of the perturbative theory as well as the non- 
perturbative theory up to the next to next to leading order 
(NNLO). We also achieve similar calculations for both the 
dispersive and the shape function models. And finally we ex¬ 
tract the coupling constant from our analysis. The last section 
summarizes our conclusions. 


II. C-PARAMETER 


The C parameter mu for electron-positron annihilation 
events is derived from the eigenvalues A/ of the linearized mo¬ 


mentum tensor 6jk. 

_ ZiPjPi/ip’i 


( 1 ) 


where p^ are the spatial components (J, k = 1 , 2 ,3) of the i-th 
particle momentum in the centre of mass frame. The sum on i 
runs over all the final state particles. If Ai are the eigenvalues 
of the matrix, we have: 


C = 3(TiT2 + T2T3+T3Ti). (2) 

The real symmetric matrix Ojk has eigenvalues At with 0 < 
/I 3 < /I 2 < Ti < 1. It describes an ellipsoid with orthogonal 
axes named minor, semi-major and major corresponding to 
the three eigenvalues. The major axis is similar but not iden¬ 
tical to ht. The C parameter varies in the range 0 < C < 1. 
C = 0 corresponds to a perfect two jet event (with massless 
jets), while C = 1 characterizes a spherical event. Planar 
events including in particular the 0(as) perturbative result, 
are distributed in the range 0 < C < 3/4. 

Eigure shows the distribution for C-parameter (0 < C < 
1) by using the Monte Carlo data (PYTHIA program) in dif¬ 
ferent center of mass energies. We see that the values are be¬ 
tween 0 and 1 for all energies. We observe a decreasing trend 
for all distributions by increasing the C parameter. The fig¬ 
ure also indicates that there is not a significant difference on 
the distributions between different energies. We conclude that 
this parameter is insensitive to the data of different energies. 


III. POWER CORRECTIONS AND DIFFERENT MODELS 

QCD calculations based on NLO perturbation theory as 
well as NNLO theory are available for many observables in 
high energy particle reactions, like the total hadronic cross 
section in e'^e~ hadrons and moments in deep inelastic 
scattering (DIS) processes and also in hadronic decay. The 
complicated nature of QCD, due to the process of gluon self 
coupling and the resulting large number of Eeynman diagrams 
in higher orders of perturbation theory, so far limited the num¬ 
ber of QCD calculations in complete NNLO 121 . We shall ex¬ 
plain the methods of calculation for NLO and NNLO theories. 
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FIG. 1: The cross section distributions for C-parameter in different 
energies. 


We also give an explanation of other ways for calculating the 
coupling constant such as the dispersive model and also the 
shape function model. These models include perturbative and 
non-perturbative regions. One of the aspects where the study 
of event shapes has taught us general lessons about QCD is 
the interplay between perturbative and nonperturbative cor¬ 
rections. The phenomenological success of renormalization 
based on models for power corrections has important conse¬ 
quences. 

In event-shape moments (y), one expects the hadronization 
corrections to be additive, such that they can be divided into a 
perturbative and a non-perturbative contribution ISl, 


All observables are chosen such that they take values between 
0 and 1. Now if we add the third sentence to above expansion, 
we will have the NNLO theory. 


1 do- _ as{ji) dA 2 dB of.Qu) 3 dC 

(Ttot dy In dy ^ In dy ^ In dry' 


C gives the NNLO correction. The perturbative contribution 
to (y") is given up to NNLO in terms of the dimensionless 
coefficients By n and Cy^n as m 


{y^)pT 




2n 


2n 


(By^n Ay nMog^—) + ( 

tLrm ^n 


2 2 2 

log^ +Ay^„(J3llog^-^ +/3ilog-^). ( 6 ) 

^y,n^ By,n and Cy^n are straightforwardly related to Byn 
and 


^y,n ~ ^y,n-> 


^y,n ~ ^y,n r) ^F^y,n ? 


C -c 

^y,n ~ ^y,n 


-CpBy n + (-C^ - K2)Ay^n- 


(7) 


</) = + <y")np- 


(3) 


Now we explicate these theories and models in the follow¬ 


ing: 


that Ay^n^ By n and Cy^n are given in Table also the constant 
K 2 is given by 191 [TOl: 


1 3 123 

K2=-^[-^CI + CfCa{—-AA^^) 


A. NLO and NNLO corrections 

The perturbative expansion of a differential distribution 
weighted of the observable y can be written for any infrared- 
safe observable for the process e^e~ hadrons. The main 
theoretical feature of jet observables is that they can be com¬ 
puted in QCD perturbation theory. As a result, jets observ¬ 
ables are finite (calculable) at the partonic level order by order 
in perturbation theory lO . If we assume the NLO theory only, 
we can compute its corresponding perturbative expansion in 
the form 

1 dcr _ asiB) . dA ^ a.Qu) 2 dB 
CFtot dy 2n dy 2n dy' 

A gives the LO result and B the NLO correction, o-tot denotes 
the total hadronic cross section calculated up to the relevant 
order. The arbitrary renormalization scale is denoted by yu. 


+ CfTrNf{-22+\6^^)1 ( 8 ) 

Here Ecm denotes the centre of mass energy squared and 
is the QCD renormalization scale. The NLO expression is 
obtained by suppressing all terms at order . The first two 
coefficients of the QCD yS-function are: 

^ llCA-4r;,A^ 


TABLE I: Contributions to C-parameter at LO, NLO and NNLO H]. 


{C') 

A 

B 

C 

n=l 

8.6379 

172.778 + 0.007 

3212.2 + 88.7 

n=2 

2.4317 

81.184 + 0.005 

2220.9 + 12.0 

n=3 

1.0792 

42.771 ± 0.003 

1296.6 + 6.7 

n=4 

0.5685 

25.816 + 0.002 

843.1 + 3.9 

n=5 

0.3272 

16.873 + 0.001 

585.0 + 3.2 
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llCl - IQCaTrNf - 6CfTrNf 
= - d - 

6 

where Ca = N color number, Cp = (N^ - l)/2N = 4/3, 
Tr = 1/2 and Nr is quark flavors. 

The perturbative coefficients By n and Cy^n) are inde¬ 
pendent of the centre-of-mass energy. They are obtained by 
integrating parton-level distributions, which were calculated 
recently to NNLO accuracy |[T|. These values are shown in 
Table U 

In order to define the non-perturbative part of the theory, 
it is necessary to explain the dispersive as well as the shape 
function model as the following. 


A 

o 



PYTHIA 

AMY 

DELPHI 

OPAL 

NLO 

•NNLO 

Dispersive Model 


100 150 

E (GeV) 


B. The dispersive model 

Non-perturbative power corrections can be related to in¬ 
frared renormalizations in the perturbative QCD expansion for 
the event-shape variables. The dispersive model for the strong 
coupling leads to a shift in the distributions ID: 

dcr dcTpt 

—(y)=^(y-ay.P) (10) 

ay ay 



where the numerical factor ay depends on the event shape that 
for C-parameter it is ay = Stt fm . Then we obtain the non- 
perturbative parameter from this model as the following: 

(/>«P = ay.P = 

- aM - (%(-) + 1 + -^) ipo a]m ( 11 ) 

Pi 4;rySo 

In the MS renormalization scheme the constant k has the 
value k = (67/18 - 7i^I6)Ca - 5I9Nf , The Milan Factor Mis 
known in two loops as M = 1.49 + 0.20 ifTTlfT^ . for number 
of flavor Nr = 3 Sit the relevant low scales ma. 

As a result for M that so-called non-inclusive Milan factor, 
we had 113: 


M = 


3.299Ca 

ySo 


+ 2 X 


-0.862Ca - 0.052Nf 

Jo 


1.575Ca-0.104A/. 

Jo 


1.49 + 0.20 


( 12 ) 


Now we can use the power corrections to calculate the per¬ 
turbative and the non-perturbative theories for the moments of 
y, we have m 


<y‘) = <j'>iVLO + (13) 

iy^) = {y^)NLO + '^{y^)NLo(ayP) + (ay.P)^, (14) 



FIG. 2: Fitting the dispersive model with data for the C-parameter 


</) = {y^)NLO + ^y^)NLo(ay.P) + 2>{y^)NLo(ay.Pf + (ay.Pf. 

(15) 

as defines the strong coupling constant and ao defines the 
nonperturbative parameter accounting for the contributions to 
the event shape below an infrared matching scale ju/ = 2. 
Therefore it is possible to obtain as and ao in this model. 
We are using the AMY data taken from TRISTAN at KEK, 
as well as DELPHI and ALEPH at CERN plus the PYTHIA 
event generator. By fitting the dispersive model with diagrams 
depicted in Eigure we observe that the model is consistent 
with the data. 

We observe that all distributions for the dispersive model 
(solid line) are in good agreement with the Monte Carlo 
(PYTHIA) and the experimental data when compared with 
NLO or NNLO prediction because the latters do not include 
the power correction. In other words the dispersive model in¬ 
cludes both perturbative and the non-perturbative part of the 
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theory. Next step is to measure the strong coupling constant 
and the non- perturbative parameter by fitting the above model 
with the data. The values for as(Mzo) and ^oOu/) up to third 
power correction are cited in Table]^ Our results are also con¬ 
sistent with those obtained by using other experiments (Table 
\nj\is extracted from CS)- 


C. The shape function model 

Korchemsky and Tafat ca describe properties of the event 
shape variables I - T and not included in NLO perturba¬ 
tion theory so called shape function, which does not depend 
on the variable nor the cms energy. This is more general than 
the dispersive model, considering both as a shift of the per¬ 
turbative prediction and as a compression of the distribution 
peak. We should mention that the shape function model in¬ 
cludes the perturbative theory as well as the non-perturbative 
part of theory. This model is a combination of both the NLO 
prediction and the power correction terms (equation |^. To 
calculate the strong coupling constant in perturbative theory, 
we are using equation]^ We also use the following expansion 
for measuring the free parameter in non-perturbative theory. 

<C‘> = <C‘)wiO + ;L. ( 16 ) 

Analogously for the second moments and the third mo¬ 
ments, we find C3: 

<C^) = {C^)nLO + '^-^{C^)nLO + (17) 

Ean 

( C ^) = {C^)nlo + ^^{C'^)nlo + ^^{C^)nlo+^- ( 18 ) 

^cm ^cm ^cm 

In order to calculate the strong coupling constant concerned 
with the perturbative theory, we perform a fitting procedure by 
including just the first part of equation On the other hand 
to find the coupling constant in the non-perturbative region we 
do a fitting procedure by taking into account both parameters 
in equation[^ We achieve a similar analysis for equations 
and The coefficient A\ is interpreted as the first moment 
and X 2 as the second moment of the shape function. These 
parameters are also known as the so called universal scales 

ca. 

Figure shows the distributions obtained from the exper¬ 
imental data as well as PYTHIA. By fitting the shape func¬ 
tion model with the corresponding distribution, we obtain the 
strong coupling constant as(Mzo) and the non-perturbative pa¬ 
rameter (Ti). Our results are indicated in Table [Tvl 


TABLE II: Measurements of the coupling constants using the disper¬ 
sive model. 


n 

1 

2 

3 

aoQui) 

0.149 + 0.003 
0.460 + 0.015 

0.117 + 0.001 
0.471 + 0.007 

0.112 + 0.003 
0.499 + 0.017 


TABLE III: Coupling constants for C-parameter in CD. 


Exp 

ALEPH 

DELPHI 

L3 

ao(Mi) 

0.123 + 0.003 
0.461 + 0.016 

0.122 + 0.004 
0.444 + 0.022 

0.116 + 0.005 
0.457 + 0.040 





EIG. 3: Eitting the C-parameter in shape function model with data. 


TABLE IV: Measurements of the coupling constants using the shape 
function model. 


n 

1 

2 

3 

/li 

0.149 + 0.003 
0.578 + 0.037 

0.143 + 0.003 
0.521 + 0.039 

0.140 + 0.0064 
0.449 + 0.033 
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These values are in good agreement with the results from 
other event shape variables carried out in our previous publi¬ 
cations Eia. Our results are also consistent with the QCD 
predictions ini. 


D. The variance 

The simple prediction for the variance of event shape vari¬ 
ables (y) on hadron level is im: 

Var(y) = {/) - {yf. (19) 



We can obtain a purely perturbative expression for the vari¬ 
ance in the dispersive model and also in the shape function 
model, up to strongly suppressed corrections 0(as(Q^)). The 
first and second orders give the identical prediction in NLO 
and NNLO moments. 

However the NLO and NNLO do not include non- pertur¬ 
bative part. 

Thus in case of the C-parameter, we have: 

Var(C) = iC^hw - iOlu)- (20) 

On the other hand, if we take into account the non- 
perturbative part of the model, we will have: 

Var{C) = {C\otai - (21) 

where the subscript indicates the perturbative as well as the 
non- perturbative parts of our calculations in both models. In 
Figures and we show the variance for dispersive as well 
as shape function model. The obtained results by fitting vari¬ 
ance with models are summarized in Table |V] and Table Ivll 
respectively. 

Both tables indicate that our values for the coupling con¬ 
stant in perturbative and in the non-perturbative regions are 
in good agreement with the QCD predictions. They are also 
consistent with those obtained from other experiments IT^ . 


FIG. 4: Fitting variance by using the dispersive model. 



FIG. 5: Fitting variance by using the shape function model. 

TABLE V: Measurements of the coupling constants by using the vari¬ 
ance of dispersive model. 


Observable 


aoiHi) 

C-parameter 

0.1039 + 0.0059 

0.5183 + 0.0211 


TABLE VI: Measurements of coupling constants by using the vari¬ 
ance of shape function model. 


Observable 


^1 

C-parameter 

0.1052 + 0.0031 

0.7164 + 0.0342 
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